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L-Ornithine cyclodeaminase (OCD) is involved in L-proline biosynthesis and catalyzes the unique 
deaminating cyclization of L-ornithine to L-proline via a A'-pyrroline^-carboxyrate (Pyr2C) inter- 
mediate. Although this pathway functions in only a few bacteria, many archaea possess OCD-like 
genes (proteins), among which only ATI 665 protein (gene) from Archaeoglobus fulgidus has been 
characterized as an NAD + -dependent L-alanine dehydrogenase (AfAlaDH). However, the physiologi- 
cal role of OCD-like proteins from archaea has been unclear. Recently, we revealed that Pyr2C reduc- 
tase, involved in trans-3-hydroxy-L-proline (T3LHyp) metabolism of bacteria, belongs to the OCD 
protein superfamily and catalyzes only the reduction of Pyr2C to L-proline (no OCD activity) [FEBS 
Open Bio (2014) 4, 240-250]. In this study, based on bioinformatics analysis, we assumed that the 
OCD-like gene from Thermococcus litoralis DSM 5473 is related to T3LHyp and/or proline metabolism 
(TILhpI). Interestingly, TILhpI showed three different enzymatic activities: AlaDH; JV-methyl-L-ala- 
nine dehydrogenase; Pyr2C reductase. Kinetic analysis suggested strongly that Pyr2C is the preferred 
substrate. In spite of their similar activity, TILhpI had a poor phylogenetic relationship to the bacte- 
rial and mammalian reductases for Pyr2C and formed a close but distinct subfamily to AfAlaDH, 
indicating convergent evolution. Introduction of several specific amino acid residues for OCD and/ 
or AfAlaDH by site-directed mutagenesis had marked effects on both AlaDH and Pyr2C reductase 
activities. The OCC_00387 gene, clustered with the TILhpI gene on the genome, encoded T3LHyp 
dehydratase, homologous to the bacterial and mammalian enzymes. To our knowledge, this is the 
first report of T3LHyp metabolism from archaea. 

© 2014 The Authors. Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies. This 
is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

L-Proline plays many important roles in the protein structure 
and cell signaling and is synthesized within cells and organisms 
from an intermediate, glutamate y-semialdehyde, commonly 
involved in biosynthesis pathways from L-glutamate and 
L-arginine: L-glutamate -» y-glutamyl phosphate -> glutamate 
y-semialdehyde (Route 1); L-arginine -> L-ornithine -> glutamate 
y-semialdehyde (Route 2) (Fig. 1A). The glutamate y-semialdehyde 
then undergoes "non-enzymatic cyclization" to give the imine, 



Abbreviations: OCD, ornithine cyclodeaminase; CRYM, u-crystallin; AlaDH, 
L-alanine dehydrogenase: Pyr2C, A'-pyrroline-2-carboxylate; T3LHyp, trans-3- 
hydroxy-L-proline; NMAlaDH, N-methyl-L-alanine dehydrogenase 
* Corresponding author. Tel./fax: +81 89 946 9848. 

E-mail address: irab@agr.ehime-u.ac.jp (S. Watanabe). 



A'-pyrroline-S-carboxyrate (Pyr5C), which is subsequently con- 
verted to L-proline by NAD(P)H-dependent Pyr5C reductase (EC 
1.5.1.2) [1]. Route 1 is the main mechanism of L-proline biosynthe- 
sis in bacteria, while eukaryotes use this route predominantly 
under stress and limited nitrogen conditions. Under high nitrogen 
input, Route 2 appears to be prominent and is mainly found in 
higher plants. 

Alternatively, L-proline is also directly synthesized from the 
amino acid L-ornithine, itself an intermediate along the above 
Route 2, by L-ornithine cyclodeaminase (OCD; EC 4.3.1.12) (Route 
3) [2]. This unique enzyme initially catalyzes NAD + -dependent oxi- 
dative deamination from the oc-amino group of L-ornithine, and 
subsequently cyclization to form A 1 -pyrroline-2-carboxyrate 
(Pyr2C) (Fig. IB). The actual cyclization step in this pathway occurs 
with the participation of enzymes, in contrast to glutamate y-sem- 
ialdehyde -> Pyr5C in Routes 1 and 2, described above. The final 
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Fig. 1. (A) Pathways of L-proline biosynthesis. (B) Proposed hydrolytic (Scheme I) and non-hydrolytic schemes (Scheme II) for the conversion of L-ornithine to L-proIine by 
OCD [9]. Step 6 corresponds to the reaction by Pyr2C reductase. Asterisks indicate putative intermediates. (C) T3LHyp pathway and the metabolic network with D-proIine. (D) 
Schematic gene clusters related to T3LHyp and/or proline metabolism of archaea (A) and bacteria (B). Putative genes in the box were purified and characterized in this study 
(see (E)). TILhpI gene was fused with OCC_00362 and OCC_00367 genes (Supplementary Methods). C-C and C-T indicate a pair of catalytic amino acid residues of proline 
racemase superfamily enzymes (see Fig. SIB). Gray putative genes are sequentially similar to other (amino acid) transporters. (E) Purification of recombinant His 6 tag 
proteins. Five micrograms each of purified protein were applied to 12% (w/v) gel (also (F)). (F) Purification of recombinant His 6 tag TILhpI mutants. 
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step by OCD is NADH-dependent reduction of Pyr2C to L-proline. 
OCD belongs to the ornithine cyclodeaminase/p-crystallin (OCD/ 
CRYM) superfamily (see Fig. 4). It was recently reported that p- 
crystallin (CRYM) from mammalians, known as an NADPH-depen- 
dent T3 thyroid hormone (triiodothyronin)-binding protein with- 
out enzymatic catalysis, functions as an NAD(P)H-dependent 
ketimine reductase (EC 1.5.1.25) [3]. Pyr2C is one of the substrates 
of this enzyme. 

Although Route 3 of L-proline biosynthesis is found in several 
bacteria, including Pseudomonas putida [2], Agrobacterium tumefac- 
iens [4], Clostridium sporogenes [5], and Sinorhizobium meliloti [6], 
many archaea possess OCD/CRYM family genes (proteins) (20- 
40% sequence identity with bacterial OCD). Among them, only 
AF1665 protein (gene) from the hyperthermophilic archaeon 
Archaeoglobus fulgidus has been characterized at the molecular 
level, and the protein functions as an L-alanine dehydrogenase 
(AlaDH; EC 1.4.1.1) but not OCD: L-alanine + NAD + + H 2 0 pyru- 
vate + NH 3 + NADH + H + [7,8]. It is proposed that the enzyme (AfA- 
laDH) is involved in L-alanine metabolism (degradation), whereas 
the catalytic efficiency value for pyruvate is much higher than that 
for L-alanine. Although OCD activity is detected in a cell-free 
extract of Methanococcus jannaschii [9], no OCD/CRYM family gene 
is found on the genome. Up to now, the physiological role of OCD- 
like protein including AfAlaDH has been unclear in archaea. 

There are two different mechanisms for OCD (Fig. IB), in which 
L-ornithine is initially oxidized to an iminium (C a = NHj contain- 
ing) intermediate (step 1 in Fig. IB). In scheme I, the iminium 
undergoes hydrolysis to 2-oxo-5-aminovalerate by water before 
cyclization (step 2), and elimination of water produces Pyr2C (step 
3). This scheme is analogous to the oxidative deamination of l- 
amino acid to 2-oxo acids by AlaDH [8], and has been preliminarily 
proposed based on biochemical experiments [10], However, the 
crystal structure and molecular dynamic analysis of OCD [2,11] 
indicate an alternative mechanism (scheme II) that no attack on 
the iminium by water occurs: intramolecular cyclic addition of C- 
5 nitrogen of the iminium to the C-2 carbon to form 2-aminopro- 
line (step 4), loss of ammonia, and reduction of the resulting imine, 
Pyr2C (step 5). 

Recently, we characterized novel bacterial pathways of trans-4- 
hydroxy-L-proline (T4LHyp) and trans-3-hydroxy-L-proline (T3LH 
yp) metabolism [12,13]. In mammalian systems, T4LHyp and 
T3LHyp occur by post-translational hydroxylation of L-proline resi- 
due in certain proteins, in particular, collagen [14]. In contrast to 
T4LHyp, the metabolic pathways of T3LHyp from bacteria and mam- 
malians may be similar [13,15] and consist of two enzymes, T3LHyp 
dehydratase (encoded by LhpH gene; EC 4.2.1.77) and NAD(P)H- 
dependent Pyr2C reductase (encoded by Lhpl gene; EC 1.5.1.1), by 
which T3LHyp is converted to L-proline via Pyr2C (Fig. 1C). Interest- 
ingly, Lhpl protein is a novel member of the OCD/CRYM superfamily, 
and catalyzes only the reduction of Pyr2C to L-proline (no OCD activ- 
ity) [11]. Furthermore, there is no sequence similarity to DpkA pro- 
tein as is known for the reductase of Pyr2C from bacteria [16], 
strongly suggesting their convergent evolution. 

Metabolism of T3LHyp (and T4LHyp) by archaea has been 
unclear until now, and it is impossible to further assess the phys- 
iological functions of OCD-like proteins in archaea based on a sim- 
ple homology search. In this study, we selected OCD-like protein 
(gene) from a hyperthermophilic archaeon, Thermococcus litoralis 
DSM 5473, as a target, and this protein functions as an NAD(P)H- 
dependent Pyr2C reductase rather than AlaDH. This gene was clus- 
tered with T3LHyp dehydratase and proline racemase genes 
involved in putative T3LHyp and/or proline metabolism. Molecular 
evolution of the OCD/CRYM family is also discussed. 



2. Results 

2.1. Candidates of target genes 

Although it is unclear whether there are archaea capable of 
metabolizing T3LHyp as a sole carbon source, a homology search 
by the Protein-BLAST program was carried out against genome 
sequences of archaea using bacterial Pyr2C reductases (Lhpl from 
Azospirillum brasilense (AbLhpI), GenBank: BA021622; Lhpl from 
Colwellia psychrerythvaea 34H (CpLhpI), GenBank: YP_268197) as 
the probe protein sequences. As described below, bacterial Pyr2C 
reductases are classified into two different subfamilies, in which 
CpLhpI and AbLhpI are contained, respectively (see Fig. 4). 
However, over 900 homologous genes (proteins) with sequence 
identities of 20-40% were found (data not shown). Similar results 
was also obtained when other OCD/CRYM family proteins includ- 
ing AfAlaDH was used instead of Lhpl protein as the probe. 

On the other hand, when bacterial T3LHyp dehydratases (LhpH 
from A. brasilense (AbLhpH), GenBank: BA021621; LhpH from 
C. psychrerythraea 34H (CpLhpH), GenBank: YP_268195) was used 
as the probe, 31 homologous genes were found. T3LHyp dehydratase 
belongs to the proline racemase superfamily, which also includes 
the archetypical proline racemase (EC 5.1.1.4) and T4LHyp epimer- 
ase (EC 5.1 .1 .8) (Fig. SI ) [ 1 5,1 7]. In contrast to bacterial Pyr2C reduc- 
tases, T3LHyp dehydratase form a single subclass in this 
superfamily. On the basis of two specific residues at the active sites, 
the members have been classified into three types: Cys-Cys type 
(proline racemase and T4LHyp epimerase); Cys-Thr type (T3LHyp 
dehydratase); Ser-Cys type (function unknown). Furthermore, two 
Cys-Cys type enzymes can be clearly distinguished in the phyloge- 
netic tree. Based on this insight, we could select only three proline 
racemase-like genes as putative T3LHyp dehydratases, although it 
has been believed that the enzyme is found only in animals and fungi 
(and bacteria) [13,15]: I litoralis DSM 5473 (OCC_00387, GenBank: 
YP_008428263 ; 45.6% identity with AbLhpH), Thermococcus sibiricus 
MM 739 (TSIB_0630, GenBank: YP_002994044; 47.1%), and Candi- 
date Korarchaeum cryptofilum OPF8 (OPF8 Kcr_0803, GenBank: 
YP_001 737233; 44.1%) (Fig. 1 ). Interestingly, these genes were clus- 
tered with OCD-like genes: T. litoralis (OCC_00362^YP_008428257 
OCC_00367, GenBank: YP_008428258); T. sibiricus (TSIB_0634, Gen- 
Bank: YP_002994048); K. ayptofilum (Kcr_0796, GenBank: 
YP_001 737226) (Fig. ID). In the cases of T. litoralis and T. sibiricus, 
another proline racemase-like gene (Cys-Cys type) and genes encod- 
ing two different types of putative heterooligomeric L-proline dehy- 
drogenase (L-ProDH) [18,19] were also contained in this gene 
cluster. This analysis strongly indicated that the gene cluster is 
related to T3LHyp and/or proline metabolism. Therefore, in this 
study, we selected GenBank: OCC_00362, GenBank: OCC_00367, 
and GenBank: OCC_00387 genes from T. litoralis as target genes 
(referred to as TILhpI and TILhpH genes, respectively). 

2.2. Preparation of recombinant His B -tag proteins 

TILhpI gene was fused with OCC_00362 and OCC_00367 genes by 
PCR (Supplementary Methods). After cloning all target genes into the 
vector pETDuet-1, the recombinant enzymes with attached His 6 -tags 
at their N-termini were expressed in Escherichia coli and purified 
with an Ni 2+ -chelating affinity column (Fig. IE). Apparent molecular 
masses of TILhpI and TILhpH, estimated by SDS-PAGE, were 40 
(36.9) and 43 (38.1) kDa (values in parentheses indicate the calcu- 
lated molecular mass of the enzyme with His 6 -tag), and those esti- 
mated by analytic gel filtration were 72 and 75 kDa, respectively 
(data not shown). Therefore, the two enzymes appear to be dimeric. 
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2.3. Characterization of TILhpI as AlaDH 

Although TILhpI shows only 44.3% sequence identity with AfA- 
laDH, we could detect both activities of the NAD + -dependent oxi- 
dative deamination of L-alanine and NADH-dependent reductive 
amination of pyruvate (in presence of 700 mM NH4): specific activ- 
ities of 0.0265 and 4.43 (unit/mg protein), respectively. NADP + was 
not used as a coenzyme. The optimum pH values of each reaction 
were 10-12 and 6-10.5, respectively (Fig. 2A). The /c cat //( m value 
for pyruvate was ~230-fold higher than that for L-alanine, caused 
by the ~120-fold lower K m value: namely, the reaction equilibrium 
favors the direction toward reductive amination (Tables 1 and 2). 
Activity was also observed with several L-amino acids and 2-oxo 
acids with a hydrophobic side chain (see below for L-proline). 
Among them, L-2-aminovalerate and 2-oxovalerate with a C 3 



aliphatic side chain were the best substrates, and the k [at /!( m values 
were 146- and 4.7-fold higher than those for L-alanine and pyru- 
vate, respectively. When pyruvate and NADH were alternatively 
incubated together with methylamine (50 mM) instead of NH4 as 
the nitrogen donor, reduction activity was detected. This result 
indicated that TILhpI functions not only as AlaDH but also N- 
methyl-L-alanine dehydrogenase (NMAlaDH; EC 1.4.1.17). Activi- 
ties of AlaDH and NMAlaDH were also detected in HPLC analysis 
using the reaction products and zymogram staining analysis 
(Figs. 2B and 3A, B).fc cat /K m values for L-alanine and pyruvate of 
TILhpI were over 1000-fold lower than those of AfAlaDH [7]. Fur- 
thermore, the K m value for NH4 of TILhpI was clearly beyond the 
physiological level (Table 2). AfAlaDH shows similar optimum pH 
of the oxidative and reductive reactions (~7.0), and no significant 
activity for hydrophobic amino acids including L-2-aminovalerate 
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Fig. 2. Enzymatic properties of TILhpI Pyr2C. (A) Effect of pH on the activity as AlaDH left panel) and Pyr2C reductase (right panel). 50 mM acetate-NaOH (pH 4.0-6.0) (square), 
50 mM potassium phosphate (pH 6.0-8.5) (circle). 50 mM Tris-HCl (pH 7.0-9.0) (triangle) for the reduction of pyruvate or Pyr2C (left axis and black symbol), and 50 mM Bis- 
tris propane (diamond) for the oxidization of L-alanine or L-proline (right axis and gray symbol), instead of 50 mM buffer under each standard assay condition. (B) Zymogram 
staining. Five micrograms each of purified protein was applied on 1 0% (w/v) gel. After electrophoresis, the gel was soaked at 50 °C (for TILhpI) or 30 °C (for AbLhpI) in staining 
solution in the presence of the indicated substrate and NAD + (each of 10 mM). Protein stained with Coomassie Brilliant Blue (CBB) was used as a loading control. (C) 
Comparison of NADH-dependent AlaDH and Pyr2C reductase activities of wild-type and mutated TILhpI proteins. The assay was performed with standard assay solution 
containing the indicated substrate (10 mM). AbLhpI and CpLhpI were assayed at 30 °C under the standard assay conditions for TILhpI. Data for AfAlaDH are from Ref. [7]. 
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Table 1 

Kinetic parameters of NAD*-dependent oxidation reaction. 



Enzymes 


Substrates 


Specific activity' (units/mg protein) 


K m (mM) 


k CM (min ') 




(min 1 mM *) 


Wild type 


L-Alanine 


0.0265 ±0.0015 


18.2 ±5.0 


2.64 ± 0.55 


0.147 


±0.012 




L-2-Aminobutyrate 


0.136 ±0.010 


6.00 ± 0.29 


6.71 ±0.16 


1.12 


±0.03 




L-2-Aminovalerate 


0.396 ± 0.030 


0.306 ± 0.009 


6.56 ± 0.08 


21.5 


±0.4 




L-Leucine 


0.248 ±0.010 


1.32 ±0.10 


10.5 ±0.4 


8.00 


±0.28 




IV-Methyl-L-alanine 


0.526 ± 0.025 


0.406 ± 0.044 


11.6 ±0.7 


28.6 


±2.1 




L-Proline 


0.384 ± 0.015 (0.0376 ± 0.0005) b 


1.12 ±0.12 (13.1 ±1.2) b 


13.8 ±1.0 (3.27±0.24) b 


12.3 


±0.4 (0.250 ± 0.006 ) b 




L-Pipecolate 


0.219 ±0.020 


0.608 ±0.115 


12.2 ±1.3 


20.3 


±1.8 


V224D/A228K 


L-Alanine 


0.0133 ±0.0010 


36.9 ±0.9 


2.12 ±0.02 


0.0575 


± 0.0008 




L-2-Aminovalerate 


0.666 ± 0.050 


2.13 ± 0.30 


7.72 ± 0.82 


3.64 


±0.12 




L-Proline 


1.51 ±0.10 


0.444 ± 0.008 


36.5 ±0.9 


82.2 


±1.1 



' Under standard assay conditions in "Experimental procedures". All substrate concentrations were 10 mM. 
b NADP*-dependent activity. 



Table 2 

Kinetic parameters of NADH-dependent reductive reaction. 



Enzymes 


Substrates 


Specific activity*' (units/mg protein) 


K m (mM) 


fc cat (min ') 


katlKn (min 'mM 1 ) 


Wild type 


Pyruvate 


4.43 ± 0.50 


9.37 ± 0.23 


315±4 


33.6 ± 0.4 




Oxalacetate 


1.78 ±0.18 


8.86 ± 0.91 


192 ±15 


21.7 ±0.5 




2-Oxobutyrate 


2.46 ± 0.25 


0.701 ± 0.020 


92.0 ±1.8 


131 ±1 




2-Oxovalerate 


2.59 ±0.21 


0.969 ± 0.093 


152 ±11 


157 ±3 




4-Methyl-2-oxovalerate 


1.02 ±0.19 


2.52 ±0.45 


129 ±18 


51.4 ±2.0 




NHJ (+2-oxovalerate) b 


2.45 ± 0.22 


1450 ±240 


267 ± 33 


0.183 ±0.007 




Pyr2C 


23.9 ±2.5 


0.944 ± 0.090 


995 ± 33 


1058 ±66 


V224D/A228K 


Pyruvate 


1.60 ±0.15 


62.0 ±12.2 


413 ±75 


6.68 ±0.11 




2-Oxovalerate 


7.00 ± 0.95 


1.82 ±0.14 


212±16 


116±0 




Pyr2C 


239 ±35 


2.77 ± 0.44 


33500 ±4900 


12900 ±300 



a Under standard assay conditions in "Experimental procedures". All substrate concentrations were 10 mM. 
b The NHJ concentrations were changed with a fixed concentration of 2-oxovalerate (10 mM). 
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Fig. 3. HPLC analysis of the reaction products. (A) Reaction product from pyruvate or 2-oxobutyrate by TILhpl (as AlaDH) in presence of NADH and NH 4 C1. Control indicates 
authentic L-alanine and L-2-aminobutyrate. (B) Reaction product from pyruvate by TILhpl (as NMAlaDH) in the presence of NADH and methylamine. Contra/ indicates 
authentic N-methyl-L-alanine. (C) Reaction product from L-proline by TILhpl (as L-ProDH) in the presence of NAD*. After the enzymatic reaction, the product was treated by 
H 2 0 2 , and then analyzed. Control indicates authentic L-glutamate and y-aminobutyrate (GABA), which are potentially yielded from Pyr5C and Pyr2C, respectively. (D) Reaction 
products from Pyr2C by TILhpl and from T3LHyp by TILhpH and TILhpl, in the presence of NADH. Control indicates authentic T3LHyp and L-proline. No peak corresponding to 
Pyr2C was observed, probably caused by no reaction with labeling reagent for amino acid analysis (data not shown). 
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(2-oxovalerate). This indicated strongly that TILhpI possesses lar- 
gely different properties from AfAlaDH, and that the physiological 
function is not the metabolism of L-alanine and pyruvate. 

2.4. Characterization of TILhpI as novel Pyr2C reductase in OCD/CRYM 
superfamily 

In contrast to AfAlaDH [7], TILhpI showed NAD + -dependent oxi- 
dization activity for L-proline. The fc cat //<m value in the presence of 
NADP + was ~50-fold lower than that in the presence of NAD + 
(Table 1): the specificity for NAD + is described above. In addition 
to L-proline, several L-proline analogs, including cis-4-hydroxy-L- 
proline (85%), trans-4-hydroxy-L-proline (70%), and L-pipecolate 
(49%), were active substrates (relative activity to L-proline in 
parentheses), whereas it is known that bacterial Lhpl proteins are 
strictly specific to L-proline [13]. These activities were also 
detected in zymogram staining analysis (Fig. 2B). When the reac- 
tion product of L-proline was further oxidized by H 2 0 2 , 4-aminobu- 
tyrate (GABA) was identified in HPLC analysis (Fig. 3C), indicating 
that TILhpI catalyzes the conversion of L-proline to Pyr2C. Gener- 
ally, L-ProDH (EC 1.5.99.8) including archaeal enzymes [18,19] is 
a flavoprotein, and the reaction product is Pyr5C (if so, L-glutamate 
is yielded by the above H 2 0 2 treatment). 

Therefore, we compared the dehydrogenation activity toward 
L-proline and other L-amino acids (Tables 1 and 2). The k CRC /K m value 
for Pyr2C (1058 min 1 mivr 1 ) was 6.7-fold higher than that for 2- 
oxovalerate, caused by the high k ca t value. Furthermore, the ratio 
of Pyr2C to L-proline in fc ca t/Xm was 86, suggesting the preference 
of the reaction equilibrium in the direction toward NADH-depen- 
dent reduction. Interestingly, the optimum pH of Pyr2C reductase 
was 4.0-6.0, clearly different from that of the reductive deamination 
of AlaDH (pH 6-10.5) (Fig. 2A). These results strongly indicated that 
Pyr2C is the preferred substrate for TILhpI. Although AbLhpI and 
CpLhpI showed only 25.3% and 31.6% sequence identity to TILhpI, 
respectively, we unexpectedly detected AlaDH activity in both 
enzymes (Fig. 3C). Commonly, these bacterial enzymes showed sig- 
nificantly high ratios of Pyr2C to pyruvate in specific activity: 176 
and 28, respectively (~5 for TILhpI). Namely, TILhpI possessed "ves- 
tigial" properties to AlaDH, probably due to the close phylogenetic 
relationship between them, as described in the main section. 

Using the reductase activity toward Pyr2C, several properties of 
TILhpI as an enzyme from "hyperthermophiles" were identified: T. 
litoralis grows between 55 and 98 °C with an optimal growth tem- 
perature of 88 °C [20]. Although the optimum temperature was 
above 75 °C, and clearly higher than that of (mesophilic) AbLhpI, 
the assay could not be performed above 75 °C due to the thermo- 
lability of NADH (Fig. S2A). Thus, TILhpI appears to be a typical 
(hyper)thermophilic enzyme. 

2.5. Phylogenetic analysis of TILhpI 

As expected from the preliminary annotation, TILhpI belongs to 
the OCD/CRYM superfamily, including the archetype OCD (subfam- 
ilies 7 and 11) [2,21], CRYM/ketimine reductase (subfamily 5) [3], 
AlaDH (subfamily 2) [7,8], L-arginine dehydrogenase (subfamily 
10) [22], L-lysine cyclodeaminase (EC 4.3.1.28; subfamily 8) [23], 
tauropine and strombine dehydrogenases (EC 1.5.1.23(22); sub- 
family 3) [24,25], L-2,3-diaminopropionate synthase (subfamily 4) 
[26], and Pyr2C reductase (from bacteria) (subfamilies 6 and 9) 
[13] (Fig. 4). TILhpI had a poor relationship to the bacterial Lhpl 
and mammalian CRYM proteins, in spite of their similar activity, 
and formed a close but distinct subfamily to many archaeal OCD- 
like proteins, including AfAlaDH (novel subfamily 1). Surprisingly, 
two OCD-like proteins from "thermophilic bacteria", Thermosedim- 
inibacter oceani (Toce_0467) and Thermovirga lienii (Tlie_1581), 
were closely related to TILhpI rather than (bacterial) AbLhpI and 



CpLhpI, suggesting the possibility of horizontal gene transfer 
between bacteria and archaea. 

Putative amino acid sequences of TILhpI contained essential 
amino acid residues for coenzyme binding (Rossmann-fold motif 
consisting of Gly-X-Gly-X 2 -[Ala/Ser], where X indicates any amino 
acid) and the catalytic tetrad for interacting (and/or binding) with 
carboxyl and amino groups of substrate ([Arg/Lys]-Lys-Arg-Asp): 
Gly 134 -Ala-Gly-Val-Gln-Ala 139 and Arg 39 -Lys 67 -Arg 110 -Asp 303 , 
respectively [2,8] (Fig. 5). Furthermore, TILhpI possessed a specific 
aspartate residue (Asp 159 ) for NAD + (H)-dependent enzymes of the 
OCD/CRYM superfamily [2,7,23-26], conforming to the coenzyme 
specificity, described above. 

2.6. Identification of catalytic amino acid residues 

It is proposed that three amino acid residues of OCD, E56-D228- 
K232 (numbering from P. putida), are related to the unique cycliza- 
tion of L-ornithine [2] (Fig. 5). In TILhpI and AfAlaDH, these amino 
acid residues correspond to L52-V224-A228 and R52- D219 - K223 , 
respectively (conserved amino acid residues are underlined). To 
obtain insight into the catalytic mechanism of OCD/CRYM super- 
family enzymes, we constructed three mutants of TILhpI, V224D/ 
A228K (DI< mutant), L52E/V224D/A228K (EDK mutant), and 
L52R/V224D/A228K (RDK mutant). They were overexpressed in 
E. coli cells as a His 6 -tagged enzyme and purified with the same 
procedures as for the wild-type enzyme (Fig. IF). Approximately 
12-fold enhancement of Pyr2C reductase activity (6.7-fold for the 
oxidization activity for L-proline), but not AlaDH activity, was 
found in the DK mutant, caused by the increasing k cat value (Tables 
1 and 2). On the other hand, when the mutations were combined 
with L52E and L52R, AlaDH activities for the 2-oxo acids were com- 
pletely eliminated, whereas Pyr2C reductase activity was the same 
as in the wild type (Fig. 2C). 

2.7. Characterization ofTlLhpH 

Potential T3LHyp dehydratase activity in TILhpH was initially 
assayed spectrophotometrically in the coupling system with TILhpI 
(in the presence of NADH). A significant decrease in absorbance at 
340 nm was observed as a result of cooperative reactions by TILhpH 
(T3LHyp -> Pyr2C) and TILhpI (Pyr2C -> L-proline): 44.6 unit/mg 
protein of specific activity with T3LHyp. K m and fc cat values for 
T3LHyp were 0.288 mM and 1710 min -1 , respectively, comparable 
with those of AbLhpH (data not shown). On the other hand, optimum 
temperature for the activity was 80-1 00 °C, which was significantly 
higher than that of the (mesophilic) AbLhpH (40-45 °C) (Fig. S2B). 
These results indicated that the TILhpH gene encodes T3LHyp dehy- 
dratase with enzymatic properties as a typical (hyper)thermophilic 
enzyme: this is the first example for archaea. 

2.8. Physiological meaning of TILhpI and TILhpH 

Although a minimal medium for the cultivation of T. litoralis is not 
avairable, we found that T. litoralis can grow in the complex medium 
[27] supplemented with L-proline, D-proline or T3LHyp as a sole car- 
bon source, instead of maltose and peptone (data not shown). On the 
other hand, Pyr2C reductase and T3LHyp dehydratase activities in 
cell-free extract prepared from T. litoralis cells grown on such as car- 
bon source were similar to those on maltose and peptone. These sug- 
gested one possibility that similar metabolic pathway of T3LHyp to 
bacteria and mammalian is operative in T. litoralis. 

3. Discussion 

The discovery of archaeal reductase for Pyr2C revealed that 
OCD/CRYM superfamily enzymes with this activity are large: 
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L-Lysine cyclodeaminase 
(Bacteria; EC 4.3.1.28) 



Streptomyces virginiae 
(BABS3673) 



Pseudomonas chlororaphis 
(Pchl06_3724) 

Azospirillum lipoferum 
(AZOLI_p40331' 



Pyr2C/Pip2C reductase 
(Bacteria; EC 1.5.1.1 (21)) 



Streptomyces pristinaespiralis 
(WP 005321563)* 
Streptomyces kanamyceticus 
(ADU56245) 



Actinoplanes friuliensis 
(AFR 18570) 



A. brasiiense 

(AbLhpl)* BurWio/ctera thailandensis 
(BTHJI1457) 

Pseudomonas aeruginosa 
(PA3862)' 



Ornithine cyclodeaminase 
(Bacteria; EC 4.3.1.12) 



Rhodobacter capsulatus 
(RCAP_rcc02113) 
Sinorhizobium meliloti 
(SM_b21494) 
Chelativorans sp. BNC1 
(Meso_3962) 



Agrohacterium tumefaciens 
(pTi-SAKURA_p190) 



Pseudomonas fluorescens 
(Pfl01_0970) 

Photobacterium profundum 
(PBPRB0030I 



Burkholderia lata 
(Bcep18194_A3938) 



Pyr2C reductase 
(Bacteria; EC 1.5.1.1) 



Ralstonia solanacearum 
(RSp0483) 

P. putida (1 U7H)* 



Shewanella violacea 
(SVI3429)" 

C. psychrerythraea 
(CpLhpl)* 



(.i-Crystallin/ketimine reductase 
(mammalian; EC 1.5.1.25) 



0.1 




L-Arginine dehydrogenase 
(Bacteria; EC 1.4.1.-) 



Catharanthus roseus 
(ABI54194)* 



Ornithine cyclodeaminase 
(Plant; EC 4.3.1.12) 



Agrobacterium tumefaciens 
(pTiBo010) 



Methanosarcina barkeri 
(Mbar_A3192)* 



Vibrio parahaemolyticus 
(VPBB 1257) 

Marinomonas mediterranea 
(Marme_3477) 



Bovine (BOS_22748) 
Human (HGNC:2418)* 
Mouse (MGM02675) 



A. fulgidus 
(AfAlaDH)* 

Archaeoglobus 
profundus (ArcpM 568} 

Thermococcus sp. 
(TES1J055)* 
Thermococcus sibiricus 
(TSIB 0634)* 

T. litoralis (TILhpl)* 



L-Alanine dehydrogenase 
(Archaea; EC 1.4.1.1) 



Herbaspirillum seropedicae 
(WP 017449912) 



ACEG dehydrogenase 
(Bacteria; EC 1.5.1.-) 



Staphylococcus aureus 
(NWMN 0061)* 

Bacillus cereus 
(WP 000933760) 



Thermosediminibacter 
oceani (Toce_0467)* 
' Thermovirga lienii 
(rile_1581 
\ Candidatus Korarchaeum 
cryptofilum (Kcr_0796) 
\ Halichondria japonica (BAD52445)" 

Suberites domuncula (CAN84630) 



Novel Pyr2C reductase 
(Archaea and bacteria; EC 1.5.1.1) 



Tauropine/strombine dehydrogenase 
(marine sponge; EC 1.5.1.23(22)) 



Fig. 4. Phylogenetic tree of OCD/CRYM protein superfamily consisting of eleven subfamilies. The number on each branch indicates the bootstrap value. Proteins with asterisks 
were used for Fig. 5A. 



subfamilies of 1, 5, 6, 7, 9, and 11 at least (Fig. 4). Furthermore, 
TILhpl (and AbLhpl and CpLhpl) possesses not only Pyr2C reductase 
activity but also (slight but significantly detectable) AlaDH and 
NMAlaDH activities. Among OCD/CRYM members, enzymatic reac- 
tions by subfamilies 2 and 10 are homologous to AlaDH, and 
enzymes of subfamilies 3 and 4 catalyze similar reactions to NMAl- 
aDH, in which taurine, glycine, and L-2,3-diaminopropionate are 
used as the nitrogen donor instead of NH4, respectively [24-26], 
Furthermore, cyclodeaminases for L-ornithine and L-lysine possess 
the same initial step of the catalytic reaction as AlaDH and NMAl- 
aDH (step 1 in Fig. IB). This suggests that the origin of L-ProDH 
(Pyr2C reductase) activity (step 6) is the same as step 1. A common 
ancestor possessed NAD(P) + -dependent dehydrogenation activity 
toward broad primary and secondary amines, and (strict) substrate 
specificity convergently evolved in bacteria, mammalians, and 
archaea. If this hypothesis is true, TILhpl acquired catalytic speci- 
ficity as Pyr2C reductase later than AbLhpl and CpLhpl, because 
the relative high activity of AlaDH and broad substrate specificity 
of L-ProDH remain (Fig. 2B and C). 

Although it is unclear whether OCD possesses (remaining) oxi- 
dative deamination activity for L-amino acid(s) except L-ornithine, 
the "potential" mechanism to produce Pyr2C (from L-ornithine) by 
TILhpl corresponds to scheme I, because of the AlaDH activity. 
Therefore, to obtain evolutional insight into Pyr2C reductase and 
OCD, we introduced several specific amino acid residues for OCD 
in TILhpl by site-directed mutagenesis. Although Asp 228 of OCD 
(numbering from P. putida) plays an important role in all steps of 
catalysis [2,11], Asp 219 conserved in AfAlaDH is not involved in 
catalysis (steps 1 and 2) [8]. Therefore, it may be reasonable that 
the V244D (and A228K) mutant of TILhpl has no significant effect 



on AlaDH activity (Tables 1 and 2, and Fig. 3C), because of the sim- 
ilar mechanism to AfAlaDH, as described above. On the other hand, 
Glu 56 of OCD favorably interacts with the 8-amino group of L-orni- 
thine. In CRYM, it is also proposed that Gly 60 at the equivalent posi- 
tion (numbering from mouse) is important for the binding of a 
large ligand, T3, because of the small side chain [28]. These insights 
explain why the substitution of Leu 52 in hydrophilic glutamate or 
arginine in TILhpl leads to a complete loss of reductive amination 
activity for hydrophobic 2-oxo acid substrates (Fig. 3C). 

As in our previous proposal [13], dehydratase activity with 
T3LHyp in the proline racemase superfamily was acquired once at 
an early evolutional stage, because TILhpH also forms a single sub- 
family together with enzymes from bacteria and mammalians 
(Fig. SI). Several (thermophilic) archaea and bacteria possess 
another proline racemase-like gene instead of the LhpH gene within 
the gene cluster containing the Lhpl gene (Fig. ID). The enzymatic 
function may be interconversion of l- and D-isomer of proline 
(and/or hydroxyproline), but not dehydration of T3LHyp: Pyr2C is 
potentially produced from D-proline by FAD-dependent D-ProDH 
(Fig. 1C). TILhpl (this study; Fig. 2B) and L-ProDH (oc 4 p 4 ) [29] show 
activities not only for Pyr2C and L-proline but also A^piperidine- 
2-carboxylate and L-pipecolate, respectively, which are involved in 
the so-called "L-pipecolate pathway" of D-lysine metabolism [16]. 
Putative L-ProDH genes are also found within the gene cluster. Fur- 
thermore, (putative) pyruvate-formate lyase (PFL) and serine/thre- 
onine dehydratase genes, producing pyruvate and/or 2- 
oxobutyrate, are often located within the flanking region of the Lhpl 
gene (Fig. ID). These indicate the possibility that the TILhpl gene 
(and the gene cluster) is also involved not only in T3LHyp but also 
in D-proline, L-proline and D-lysine metabolism, and several 2-oxo 
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Fig. 5. (A) Partial multiple sequence alignment of deduced amino acid sequences of TILhpI. Binding sites for a carboxyl group of substrate are shaded in violet, green, yellow 
and light-green. Aspartate, lysine and glutamate residues, related to the cyclization by OCD, are shaded in blue, orange, and red, respectively. A specific aspartate residue for 
NAD*(H)-dependent enzymes, potentially interacting with 2'- and 3'-hydroxyl groups in the ribosyl moiety of NAD + (H), are shaded in pink. Amino acid residues Coenzyme- 
binding motif of Rossmann-fold are shown as white letters in black boxes. Gray-shaded letters indicate highly conserved amino acid residues. (B) Schematic diagram showing 
the interactions of [.-ornithine (left panel) and L-alanine (right panel), and nearby residues in OCD and AlaDH, respectively. Color of residues correspond to Fig. 5A. (For 
interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 



acids. Although T. litoralis can grow on T3LHyp as a sole carbon 
source, Pyr2C reductase (and T3LHyp dehydratase) activity is not 
significantly induced in cell-free extract, in contrast of bacteria 
[13], conforming to the putative functions in multiple metabolism 
described above. Development of gene disruption of T. litoralis 
would be useful for understanding more detailed physiological 
role(s) of TILhpI gene. 

4. Experimental procedures 

4.1. General procedures 

Basic recombinant DNA techniques were performed as described 
by Sambrook et al. [30]. PCR was carried out using a GeneAmp PCR 
System 2700 (Applied Biosystems) for 30 cycles in a 50 |il reaction 
mixture containing 1 U of KOD FX DNA polymerase (TOYOBO), 
appropriate primers ( 1 5 pmol) and template DNA under the follow- 
ing conditions: denaturation at 98 °C for 10 s, annealing at 50 °C for 
30 s and extension at 68 °C for the time periods calculated at an 
extension rate of 1 kbpmin -1 . DNA sequencing was carried out 
using the BigDye Cycle Sequencing Kit ver.3.1 (Applied Biosystems) 
and appropriate primers with the Genetic Analyzer 3130 (Applied 
Biosystems). Protein concentrations were determined by the 
method of Lowry et al. [31 ] with bovine serum albumin as the stan- 
dard. SDS-PAGE was performed as described by Laemmli [32]. 

4.2. Strain and growth conditions 

T. litoralis DSM5473 were grown under anaerobic conditions at 
80 °C in a complex medium, as described previously [27], in the 



absence of sulfur. If necessary, maltose, peptone and yeast extract 
were omitted, and L-proline, D-proline or T3LHyp (30 mM) was 
added alternatively. 

4.3. Plasmid construction for expression of recombinant proteins 

Primer sequences used in this study are shown in Table SI. In 
this report, the prefixes Tl (T. litoralis), Af (A. fulgidus), Ab (A. brasi- 
lense) and Cp (C. psychrerythraea) have been added to gene symbols 
or protein designations when required for clarity. TILhpI (fusion of 
OCC_00362 and OCC_00367; see Supplementary Methods) and 
TILhpH genes (OCC_00387) were amplified by PCR using primers 
containing appropriate restriction enzyme sites at the 5' and 3' 
ends and genome DNA of T. litoralis as a template. Each amplified 
DNA fragment was introduced into BamHl-Pstl sites in pETDuet- 
1 (Novagen), a plasmid vector for conferring N-terminal His 6 -tag 
on expressed proteins, to obtain pET/TlLhpIwT and pET/TlLhpH. 

4.4. Expression and purification of His 6 -tagged recombinant proteins 

E. coli BL21(DE3)-RIL (Novagen) harboring the expression plas- 
mid for His 6 -tagged proteins was grown at 37 °C to a turbidity of 
0.6 at 600 nm in Super broth medium (pH 7.0, 12 g tryptone, 
24 g yeast extract, 5 ml glycerol, 3.81 g KH 2 P0 4 , and 12.5 g 
K 2 HP0 4 per liter) containing 50 mg/liter ampicillin. After the addi- 
tion of 1 mM isopropyl-|3-D-thiogalactopyranoside (IPTG), the cul- 
ture was further grown for 6 h at 37 °C to induce the expression 
of His 6 -tagged protein. Cells were harvested and resuspended in 
Buffer A (50 mM sodium phosphate buffer (pH 8.0) containing 
300 mM NaCl and lOmM imidazole). The cells were then 
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disrupted by sonication for 20 min at appropriate intervals on ice 
using the Ultra Sonic Disruptor Model UR-200P (TOMY SEIKO Co., 
Ltd., Tokyo, Japan), and the solution was centrifuged at 
108,000xg for 20 min at 4 °C. The supernatant was loaded onto a 
Ni-NTA Superflow column (Qiagen) equilibrated with Buffer A 
linked to the BioAssist eZ system (TOSOH). The column was 
washed with Buffer B (50 mM sodium phosphate buffer (pH 8.0) 
containing 300 mM NaCl, 10% (v/v) glycerol, and 50 mM imidaz- 
ole). The enzymes were then eluted with Buffer C (pH 8.0, Buffer 
B containing 250 mM imidazole instead of 50 mM imidazole), con- 
centrated by ultrafiltration with Centriplus YM-30 (Millipore), dia- 
lyzed against 50 mM Tris-HCl buffer (pH 8.0) containing 50% (v/v) 
glycerol, and stored at -35 °C until use. If necessaiy, AbLhpI and 
CpLhpl proteins were prepared, as described previously [13]. 

The native molecular mass of recombinant proteins was esti- 
mated by gel filtration, which was carried out using HPLC with a 
Multi-Station LC-8020 model II system (TOSOH) at a flow rate of 
1 mL min -1 . The purified enzyme (~10 mg mlr 1 ) was loaded onto 
a TSKgel G3000SWXL column (TOSOH) equilibrated with 50 mM 
Tris-HCl buffer (pH 8.0). A high molecular weight gel filtration 
calibration kit (GE Healthcare) was used as a molecular marker. 

4.5. Enzyme assay 

All enzyme assays were performed at 50 °C unless otherwise 
indicated. 

Pyr2C reductase was assayed routinely in the direction of Pyr2C 
reduction by measuring the oxidization of NADH at 340 nm, using 
a Shimadzu UV-1800 spectrophotometer (Shimadzu GLC Ltd., 
Tokyo, Japan). The standard assay mixture contained 10 mM Pyr2C 
in 50 mM potassium phosphate (pH 6.0) buffer. The reactions were 
started by the addition of 100 ul of a 1.5 mM NADH solution to a 
final volume of 1 ml. To assay the reverse reaction, the reaction 
mixture consisted of 50 mM Bis-tris propane (pH 10.5) and 
10 mM L-proline. The reaction was started by the addition of 
15 mM NAD + (100 ul) with a final reaction volume of 1 ml. One 
unit of enzyme activity refers to 1 umol NADH produced/min. 
Pyr2C was enzymatically synthesized from T3LHyp, described pre- 
viously [13]. 

AlaDH activity was measured as oxidative deamination of L-ala- 
nine and as reductive amination of pyruvate. For the oxidative 
deamination assay, the standard assay mixture contained 50 mM 
Bis-tris propane (pH 10.5), 10 mM L-alanine, and 1.5 mM NAD + . 
The reductive amination assay mixture consisted of 50 mM Tris- 
HCl (pH 8.0), 700 mM NH 4 C1, 10 mM pyruvate, and 0.15 mM 
NADH. For the assay of NMAlaDH activity, 50 mM methylamine 
was used instead of 700 mM NH 4 C1. 

T3LHyp dehydratase was assayed spectrophotometrically in the 
coupling system with TILhpI. The reaction mixture consisted of 
50 mM Tris-HCl (pH 8.0), 0.15 mM NADH, and 10|ig purified 
TILhpI. The reaction was started by the addition of 100 mM T3LHyp 
(100 with a final reaction volume of 1 ml. The enzyme was 
alternatively assayed by the colorimetric method based on the 
reaction of 2-aminobenzaldehyde with Pyr2C, which yields a yel- 
low reaction product [15]. This method was used to determine 
the optimum temperature for activity. 

4.6. Reaction product analysis 

Purified TILhpI (10ug) was added to the following four mix- 
tures (1 ml): for AlaDH activity, 50 mM Tris-HCl (pH 8.0) buffer 
containing 10 mM pyruvate (or 2-oxobutyrate), 500 mM NH 4 C1, 
and 10 mM NADH; for NMAlaDH activity, 50 mM Tris-HCl (pH 
8.0) buffer containing 10 mM pyruvate, 100 mM methylamine, 
and 10 mM NADH; for L-ProDH activity; 50 mM Bis-tris propane 
(pH 10.5) buffer containing 10 mM L-proline and 10 mM NAD + ; 



for T3LHyp dehydratase activity of TILhpH; 50 mM Tris-HCl (pH 
8.0) buffer containing 10 mM T3LHyp, 10 mM NADH, and TILhpH 
(10 ug). After incubation at 50 °C overnight, each enzyme product 
was then analyzed by a Hitachi L-8900 amino acid analyzer (Tokyo, 
Japan), using ion exchange chromatography followed by post-col- 
umn derivatization with ninhydrin. In the case of L-ProDH activity, 
the product was further oxidized by the addition of 30% H 2 0 2 
(20 ul) to the mixture before HPLC analysis [33]. 

4.7. Zymogram staining analysis 

Purified TILhpI was separated at 4 °C on non-denaturing PAGE 
with 10% (w/v) gel, which was performed by omitting SDS and 2- 
mercaptoethanol from the solution used in SDS-PAGE. The gels 
were then soaked in 10 ml staining solution consisting of 50 mM 
Bis-tris propane (pH 10), 0.25 mM nitroblue tetrazolium (NBT), 
0.06 mM phenazine methosulfate (PMS), 10 mM substrate, and 
10 mM NAD + at 50 °C for 5 min in the dark. Dehydrogenase activity 
appeared as a dark blue band. 

4.8. Site-directed mutagenesis 

The mutation in TILhpI was introduced by sequential steps of 
PCR [34] using sense and antisense primers (Table SI) and pET/ 
TlLhpIwr (for V224D/A228K (DK) mutant) or pET/TlLhpI DK (for 
L52E/V224D/A228K (EDK) and L52R/V224D/A228K (RDK) 
mutants) as a template. Briefly, the codons used for single mutants 
were as follows: L52E, CTG -> AGG; L52R, CTG -> GAG; V224D, 
GTG -» GAC; A228K, GCC -> AAG. The coding region of the mutated 
genes was confirmed by subsequent sequencing in both directions. 
Mutant proteins were expressed and purified by the same proce- 
dures as for the wild-type enzyme. 

4.9. Amino acid sequence alignment and phylogenetic analysis 

Protein sequences were analyzed using the Protein-BLAST and 
Clustal W program distributed by DDBJ (DNA Data Bank of Japan) 
(www.ddbj.nig.ac.jp). Multiple sequence alignment was performed 
with the following default parameters: Protein Weight Matrix, Gon- 
net; GAP OPEN, 10; GAP EXTENSION, 0.20; GAP DISTANCES, 5; NO 
END GAPS, no; ITERATION, none; NUMITER, 1 and CLUSTERING, 
NJ. The phylogenetic tree was constructed by the neighbor-joining 
method [35]. Bootstrap resampling was performed 1000 times. 
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